NACA TN 2763 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


GUST-RESPONSE ANALYSIS OF AN AIRPLANE INCLUDING 

WING BENDING FLEXIBILITY 

By John C. Houbolt and Eldon E. Kordes 

Langley Aeronautical Laboratory 
Langley Field, Va. 


TECHNICAL NOTE 2763 



Washington 
August 1952 


IE 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2763 

GUST-RESPONSE ANALYSIS OF AN AIRPLANE INCLUDING 
WING BENDING FLEXIBILITY 1 
By John C. Houbolt and Eldon E. Kordes 

SUMMARY 


An analysis is made of the gust response (including bending moment) 
of an airplane having the degrees of freedom of vertical motion and wing 
bending flexibility and basic parameters are established. A convenient, 
but accurate, numerical solution of the response equations is developed 
which is very well suited for making trend studies. An example treated 
shows results which are in very good agreement with the results obtained 
by a more precise but more lengthy method. 

The method of determining a gust causing a known response is indi- 
cated and a procedure is given for determining the response of an air- 
plane directly from the known response of another airplane by eliminating 
the common gust condition. 


INTRODUCTION 


In the design of aircraft the condition of gust encounter has become 
critical in more and more instances, mainly because of the ever-increasing 
flight speeds. Aircraft designers have therefore placed greater emphasis 
on obtaining rational methods for accurately predicting the stresses that 
develop. As a result, the number of papers dealing with the prediction 
of stresses in an aircraft traversing a vertical gust has significantly 
increased. (See, for example, refs. 1 to 9* ) Many of the papers have 
treated the airplane as a rigid body and in so doing have dealt with 
either the degree of freedom of vertical motion alone (refs. 5 to 8) or 
with the degrees of freedom of vertical motion and pitch (refs. 7 and 9). 


■^This paper is a revision and extension of a paper entitled n The 
Determination of the Response Due to Gusts of One Airplane From the Known 
Response of Another Airplane” published as TN No. Structures 40, British 
R.A.E., June 19^9* which was completed by Mr. Houbolt during a temporary 
tour of duty with the Royal Aircraft Establishment. Since the present 
paper is complete in itself, no further reference to the earlier paper 
is necessary. 
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Of greater concern in the consideration of gust penetration, however, is 

the influence that wing flexibility has on structural response. This 

concern has two main aspects: (l) that including wing flexibility may 

lead to the calculation of higher stresses than would be obtained by 
rigid-body treatment of the problem and (2) that wing flexibility may 
introduce some error when the airplane is used as an instrument for 
measuring gust intensity. Thus, there are also many recent papers which 

treat the airplane as an elastic body. In most of these papers the 

approach used involves the development of the structural response in 
terms of the natural modes of vibration of the airplane (refs. 2 to 4). 
Others used a more unusual approach, as for example, reference 1, which 
deals with the simultaneous treatment of the conditions of equilibrium 
between aerodynamic forces and structural deformation at a number of 
points along the wing span. Whatever the approach, however, the main 
disadvantage of these elastic-body analyses is that they are not very 
well suited for making trend studies without excessive computation time. 

In the present paper, the case of the gust penetration of an air- 
plane having the degrees of freedom of vertical motion and wing bending 
is considered. Wing bending was chosen because designers have expressed 
greater concern about the influence of this flexibility on gust response 
than they have about other types of flexibility. The paper has the 
objective of trying to establish some of the basic parameters that are 
involved when wing bending flexibility is included and of developing a 
method of solution which is fairly well suited for trend studies without 
excessive computation time. Such a procedure would be useful in evalu- 
ation studies which are intended, for example, to evaluate the effect of 
such factors as forward speed, spanwise mass distribution, gust length, 
and gust shape. 

The equations for response (including accelerations, displacements, 
and bending moments) are derived and the basic parameters outlined. An 
easy numerical solution for the response which is readily handled either 
by manual or machine methods is then given. The inverse of the response 
problem is considered briefly; that is, the method of determining the 
gust causing a given response is indicated. Finally, on this basis, a 
procedure is outlined whereby the response of one airplane may be found 
directly from the known response of another airplane without the 
necessity of establishing the gust causing the known response. 


SYMBOLS 


Tl 


slope of lift curve 


deflection coefficient for nth mode, function of time alone 
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A 

b 

c 

c o 

E 

f(s) 

F 


H 

I 

K J 

Iry 

L g 

m 

M J 

M c 0 

M n 

An 

P 

s, a 
S 


aspect ratio of wing 

span of wing 

chord of wing 

chord of wing midspan 

Young’s modulus of elasticity 

nondimensional gust force, / - cr)dcr 

J 0 u 

external applied load per unit span 
acceleration due to gravity 
distance to gust peak, chords 
bending moment of inertia 

nondimensional bending-moment factor = Kj — pUVM c ^ 

aerodynamic lift per unit span of wing due to vertical motion 
of the airplane 

aerodynamic lift per unit span of wing due to gust 
mass per unit span of wing 

net incremental bending moment at wing station j 

moment of wing area about spanwise station under consideration 

generalized mass of nth mode 

incremental number of g acceleration 

load intensity per unit spanwise length 

distance traveled, 2Z. t, half-chords 

c o 

wing area 
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time, zero at beginning of gust penetration 
vertical velocity of gust 
maximum vertical velocity of gust 
forward velocity of flight 
total weight of airplane 

distance along wing measured from airplane center line 

deflection of elastic axis of wing, positive upward 

deflection of elastic axis in nth mode, given in terms of a 
unit tip deflection 

response coefficient based on a n , ^ — a n 

second derivative of Zq with respect to s 

second derivative of with respect to s 

bending-moment response factor, ratio of bending moment 

obtained for airplane considered flexible to bending moment 
obtained for airplane considered rigid 

distance interval, half-chords 


reduced- frequency parameter, 


^o 

2V 


nondimens ional relative-density parameter. 


8M n 

apc Q S 


nondimens ional bending-moment parameter, - — 

apc o M C0 

mass density of air 

function which denotes growth of lift on rigid wing entering 
a sharp-edge gust (Kussner function) 
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natural circular frequency of vibration of nth mode 

1 _ <x> function which denotes growth of lift on an airfoil following 
a sudden change in angle of attack (Wagner function) 

e = 1 - ® 

Subscripts: 

j spanwise station 

n natural modes of vibration 

m number of distance intervals traveled 

Notation: 


column matrix 



square matrix 


Dots are used to denote derivatives with respect to time; primes 
denote derivatives with respect to s or 0 . 


ANALYSIS 

Equations for Structural Response 


Equations of motion .- Consider an airplane flying horizontally into 
vertical gusts, and suppose that it is desired to include wing bending 
flexibility in determining the stresses induced by these gusts. The 
problem is actually one of determining the response of an elastic wing 
subject to dynamic forces. For dynamic forces of intensity F per unit 
length, the differential equation for wing bending is, if structural 
damping is neglected, 


El ^ = -mw + F (1) 

£>y 2 dy 2 


where w is the deflection of the elastic axis referred to a fixed 
reference plane. The task of determining the deflection that results 
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from the applied forces F may be handled conveniently by expressing 
the deflection in terms of the natural free- free vibration modes of the 
wing. With regard to the flight of an airplane through gusts, exami- 
mation of a number of acceleration and strain records that have been 
taken in normal flight with several different aircraft shows that the 
response to gusts is composed primarily of a rigid-body vertical trans- 
lation and fundamental-bending-mode excitation of the wing. Thus, the 
assumption is made in the present analysis that the response may be given 
with fair accuracy by considering only these two degrees of freedom. 

This assumption is probably invalid when the airplane is flying near the 
flutter speed, for then a large amount of coupled bending-torsion dis- 
placement may occur. (See ref. 3 *) 

The wing deflection is thus assumed to be given by the equation 


w = a Q + a lWl 


( 2 ) 


where w^ is the deflection given in terms of a unit tip deflection 
along the elastic axis of the wing for the fundamental mode, and a Q 
and a]_ are functions of time alone. In this form aQ denotes the 

free-body vertical displacement of the airplane (in this case the dis- 
placement of the nodal points) and a-|_ is the part of the wing-tip 

deflection which is associated with the fundamental mode, as illustrated 
in the following sketch: 



The use of symmetrical modes implies that only the symmetrical gust is 
to be considered hereinafter. 

Substitution of equation (2) into equation (l) yields 


El 


dy' 


d 2 w-L 

dy 2 


= -m 


( S 0 + Vl) 


a 


( 3 ) 
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From the following relation which expresses the condition for natural 
fundamental-mode vibration 




O 

cc^mwj 


equation ( 3 ) may be written 


a^cD^'Tnw-^ = -m(aQ + ^l w l) + F 


(4) 


where o>j_ is a natural circular frequency of vibration of the funda- 
mental mode. If this equation is integrated over the wing span and use 
is made of the following known orthogonality condition of the free-body 
and fundamental modes: 



mw-|_ dy 


0 


the following equation results: 


rW 2 

M o a o = / F d y 

J -b/2 


(5) 


( 6 ) 


where Mq is the airplane mass. Now, if equation (4) is first multi- 
plied through by w^ and then integrated over the wing span and use is 
made of equation (5), the following equation is obtained: 


2 r h/2 

M^a^ + = / Fw-^ dy ( 7 ) 

-b/2 

where is the generalized mass for the fundamental mode, that is, 

f b / 2 ? 

M-|_ = / mw-^dy. Equations (6) and (7) represent, respectively, 

J -b/2 

the equations for free-body motion and fundamental wing bending and can 
be solved if the forces F are known. 
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For the present case of the airplane flying through a gust, the 
force F is composed of two parts: a part designated by Ly due to 

the vertical motion of the airplane (including both rigid-body and 
bending displacements) and a part Lg resulting directly from the gust 
(this latter part is the gust force which would develop on the wing 
considered rigid and restrained against vertical motion). These two 
parts are defined (see refs. 1 and 3) in the equation of F as follows: 

"t o"t 

F = Ly + Lg =-| pcvj" w[l - <t(t - T^dT + | pcV J 0 ut(t - T)dT (8) 


where 1 - <t(t) is a function (commonly referred to as the Wagner 
function) which denotes the growth of lift on a wing following a sudden 
change in angle of attack and for two-dimensional incompressible flow 
is given by the approximation 


& 


-o.09^t -0.6jt 

<I>(t) = 1 - 0.l65e - 0.335e 

-* A=<» 


(9) 


and \lr(t) is a function (commonly referred to as the Kussner function) 
which denotes the growth of lift on a rigid wing penetrating a sharp- 
edge gust and for two-dimensional incompressible flow is given by the 
app r oxima t i on 



-0.26^-t 

0. 5e 


0.5e 



( 10 ) 


An additional term which involves the apparent air mass should be included 
in equation (8); this mass term is inertial in character and may be 
included with the structural mass (see ref. 1) although it is usually 
small in comparison. The lift-curve slope a may be chosen so as to 
include approximate over-all corrections for aspect ratio and compressi- 
bility effects. 

If w as given by equation (2) is substituted into equation (8) 
and the resulting equation for F is substituted into equations (6) 
and (7), the following two equations are obtained for the case of a 
uniform spanwise gust: 


2M 0 

apVS 



Zl 

S 


apjfl - ®(t - T)]d' 



ui|r(t “ t H t (ID 


2E 
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2M-) .. co, 2 2M 1 

— — a, + — a n 

apVS x apVS 1 



7 *l) C 1 ' t(t - t D 


dT + 



t)1t 


( 12 ) 


where (because of mode symmetry) 




ob/2 

S = 2 / c dy 

J 0 


nb/2 

s i = 2 / cw i dy r 


-b/2 


S 2 = 2 / cw^ dy 


(13) 


Equations (ll) and (12) may be put in convenient nondimensional 
form by introducing the notation 


s 



or 



(14) 



(15) 


where c Q is the midspan chord of the wing and U is the maximum 
vertical velocity of the gust. With this notation, equations (ll) 
and (12) may be written 
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^1 Z 1 


d -^ 2 z ± 


= -2 


, r i z o 


TpZ 


[l - $(s - cruder 



t(s - a)da 


(17) 


where 


8Mp 

" apc Q S 


8M-, 

LL = = — 

1 apc G S 


X = 


2V 





(18) 


and a prime denotes a derivative with respect to a. Equations (l6) 
and (17) are the basic response equations in the present analysis. The 
five parameters appearing in these equations and given by equations (18) 
depend upon the forward velocity, air density, lift-curve slope, and the 
airplane physical characteristics: the wing plan form, wing bending 

stiffness, and wing mass distribution. Experience has shown that vari- 
ations in the physical characteristics cause significant variations in 
the first three of the five parameters, while the last two vary only to 
a minor extent. The first three are therefore the most basic param- 
eters; Pq is a relative-density factor, frequently referred to as a 

mass parameter, and is associated with vertical free-body motion of the 
airplane; p^_, similar to Pq, is the mass parameter associated with the 

fundamental mode; and X by its nature may be interpreted as a reduced- 
frequency parameter similar to that used in flutter analysis. 
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It is significant to note that, if any one of the three quant- 
ities z Q , z-j_, and u appearing in equations (l6) and (17) is specified 

or known, the other two may be determined. Thus, if the gust is known, 
the response may be determined, or conversely, if either z Q or z 1 is 

known, the gust may be determined. A useful equation relating Zq and 
z may be found by combining equations (16) and (17) so as to eliminate 
the integral dealing with the gust. The result is the equation 

7^( z l" + x,2z i) + ~ r ijj' z l"[' 1 - “ a )] da = dQ z o" ^ 1 9) 

which is used subsequently. 

It may also be of interest to note that ^ 0 z o" in effect d - efines a 
frequently used acceleration ratio. From equations (12) and (11), the 
rigid-body component of the vertical acceleration may be written 



z 0 


It 


or, when expressed in terms of the incremental number of g f s, 


An 


_ a 0 _ 4VU _ t, 

S c 0 g 0 


An acceleration factor An s based on quasi-steady flow and peak gust 
velocity is now introduced according to the definition 

An s w - | psv 2 u. 
is thus found to be 

An s 



^o z o" 


The ratio 
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Where the gust shape is represented analytically and the unsteady- 
lift functions are taken in the form given by equations (9) and (10), 
solution of the response equations may he made by the Laplace transform 
method, but such a solution is more laborious than desired. Therefore, 
a numerical procedure which permits a rather rapid solution of the equa- 
tions has been devised and is presented in a subsequent section. It may 
be well to mention, however, that the response equations are suitable 
for solution by some of the analog computing machines. 

Bending stresses .- The bending moment and, hence, the bending 
stresses that develop in the wing due to the gust may be found as 
follows: The right-hand side of equation (1) defines the loading on 
the wing; suppose this loading is noted by p, then 


P 


-mw + F 


By use of equations (2) and (8), and the notation of equations (l4) 
and (15), this equation becomes 

p = -m ^U(z 0 " + z 1 M w 1 ) - apcVuT (z 0 " + z 1 "w 1 ) fl - $>(s - a)l da + 
c o ^ 0 L J 



u ! ^(s - a)da 


If the moment of this loading is taken about a given wing station, 

say y the following equation for incremental bending moment at that 
J 

station would result: 

rW 2 

M j =J p ( y - y j) dy 

y j 


- +Mm i zi ") ■ ap ™/T( 


M c z 0 + M c n z l 


")[l 



da + 


a 

2 


PVMcq 


0 


u 1 s 


a)da 


,s 


( 20 ) 


NACA TN 2763 


13 


where the M's bearing double subscripts are first moments defined as 
follows : 



an< ^ being considered. Division through of eq.ua- 

tion (20) by the quantity | pVUM C() gives the following equation which 
is considered to define a bending-moment factor Kj at wing station y 


K i 


J__ 

1 pVUM c 

2 c 0 


. . J z ■ 

apcoM Co ^ 0 






da + 



a)da 


( 22 ) 


The factor || pVUMc^ 


may be regarded as the maximum aerodynamic bending 


moment that would he developed by the gust under conditions of quasi- 
steady flow and with the wing considered rigid and restrained against 
vertical motion at the root. The bending-moment factor Kj may thus 
be seen to be the ratio of the actual dynamic bending moment that occurs 
to this quasi-steady bending moment and therefore may be regarded as a 
response or an alleviation factor. 
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A more convenient form for the bending-moment factor may be 
obtained by solving equations (l6) and (17) simultaneously for the 

quantities f z Q " [I - *(s - a)] da and f zp" [l - <t(s - a)]dc and 

Jo J 0 

substituting these values into equation (22). With these operations the 
following equation results: 


M 1 

K i = “ 

J | pVUM c 


„ 0 . n h 0 ” ♦ H - nVi" * -l^l 


■ r l ' r 2 


where 


r l - r 2 


r l - r 2 


(23) 


M C] 
r 3 = M c 


8 V 


^0 apc^ 


8m, 


l m 


= 


1 apCoM c 




(24) 


It is seen that, when bending moments are being determined, three addi- 
tional basic parameters (eqs. (24)) appear. The similarity of 
and q 1 to p Q and Up is to be noted; first moments of masses and 
areas are involved rather than masses and areas. 

Reduction to rigid case .- It may be of interest to show the reduc- 
tion of the response equation to the case of the airplane considere as 
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a rigid body. Thus, if z ^ is equated to zero in equation (16), the 
following equation for rigid-body response is obtained: 



If z^" is set equal to zero in 

equation (25), the following equation 
for the rigid-body case is obtained 

Kj = (HQ - 



equation (22) and use is made of 
for the bending-moment parameter 

7o) z o" ( 2 6) 


where Zq" is the acceleration of the airplane considered as a rigid 
body. 


Matrix Solution of Response Equations 

In this section a rather simple numerical solution of the response 
equations (l6) and (17) is presented. The procedure is readily adapted 
to either manual or punch-card machine calculations. 


The derivation proceeds on the basis that the response due to a 
given gust is to be determined. The airplane, just before gust pene- 
tration, is considered to be in level flight and hence has the initial 
conditions that the vertical displacement and vertical velocity are both 
zero. These conditions mean that Zq, z^, Zq' , and z-^* are all zero 

at s = 0. The gust force can be shown to start from zero and, there- 
fore, the additional initial conditions can be established that Zq" 

and z-^ are also zero at s = 0. By the numerical procedure, solution 

for the response at successive values of s of increment e will be 
made and, for the case being considered, it is found advantageous to 
solve directly for the accelerations rather than the displacements. 

In order to make the presentation more compact, the following 
notation is introduced: 


a 


s 0 


P = 




(27a) 


0 = 1 -$ 


1 6 
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and 


f(s) = J I|r(s - a)da 


(27b) 


With this notation, equation (l6) would appear simply as 


P 0 a = -2 / + r-^P>\9(s - cr)da + f(s) 

Jq v y 


( 28 ) 


In accordance with numerical-evaluation procedures, the interval 
between 0 and s is divided into a number of equal stations of 
interval e so that s = m€. The product of (a + r ]_P) and 0(s - a) 

is assumed formed at each station and, with the use of the trapezoidal 
method for determining areas, the unsteady-lift integral in equation (28) 
may be written in terms of values of a and (3 at successive stations 

as follows, where the mth station corresponds to the value s: 

( a + r x P)^(s - tf)da = 6 ^ 0 m-l a l + 0 m-2 a 2 + • • • + 0 l a m-l + \ + 

er l(^m-l^l + ^m-2^2 + + ^l^m-1 + \ ^O^m) 

(29) 

in which 0 q, 9 ±, . . . are, respectively, the values of the 
1-0 function at s = 0, s = e, . . . (o,q and do not appear 

because of the initial conditions). With this equation, equation (28) 
may be written at various values of s or at successive values of m; 
the result, for example, for m = 1 is 

p 0 ai = -^o a l “ ^ r l 9 oh + f l 

and for m = 2, 

0 

[i Q a 2 = -€(20^ + d(f- 2 ) ~ 6r l( 20 lPl + 0 O P 2) + f 2 


CO 

M 


where and f2 are the values of the gust-force integral at s = e and s = 2e. The 

equations thus formed may he combined in the following matrix equation: 


n 0 + 0o € 



a l 


r l e 0 € 




Pi 


f l 

2 

H 0 + 0 0 € 


a 2 



ri0 o € 



?2 


f 2 

2^2^ 

20 x e + 0 q 6 

• . 


a 3 

+ 

2r 1 6 2 e 

^ r l®l e 

r i d cf 


p 3 

= 

f 3 

20 

m-1 

20 m-2 e * * * ^O +0 O e 


a m 


2r-,0 rn ,€ 
1 m-1 

2r l 0 m-2 e 

• • • r l 0 O € 


Pm 


f m 


(30a) 


which may be abbreviated 


HM + HM - l fi 


(30b) 


The simplicity of the matrices A and B, and all square matrices to follow, is to be noted; 
the matrices are triangular and all elements in one column are merely the elements in the 
previous column moved down one row. Thus, only the elements in the first columns have to be 
known to define completely the matrices. 

Now instead of considering directly the second response equation, equation (17) > it is 
expedient to consider equation (19) which is repeated here for convenient reference 


|l(zi" + (Aj ♦ 2^£ - r ^ ^”8(8 - c)do 


■ Vo 
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According to the derivation presented in the appendix, the value of z ]_ at s me may be 
approximated in terms of the past-history value of by the following equation: 


z 1 = e 2 
J -m 


(m - 1)P]_ + • • • + 2 Pm-2 + ^m-l + \ ^m 


( 31 ) 


where P]_, 0 2 > • • • are the values of z l" at s " e > s - 2e , . . . . If th ls equation is 

used to replace Z] _ in equation (19) and the unsteady-lift integral is manipulated similarly 

to the integral in equation (28), equations are obtained for successive values of m which 
involve only the unknowns a and p. The results may be combined in the following matrix 

equation: 



— 




^( 1 + 4 £ ) + @- r i ) e ° 6 




a l 

^xV + 2 g- ri ) V 


p 2 


a 2 

2 111 x 2 6 2 + 2g- r 1 )e 2 € 


p 3 


a 3 

(m-lg ^ +2 g r l)Vl e 

i.-2)Sl A 2 .2g- ri )vi* 7( 1 * ! ¥)*(7' r ‘) Srj ' 





( 32 a) 


DOM = ^ol a 


which may be written 


( 32 b) 
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The square matrix C is seen to be similar to the other square 
matrices in that it is triangular with all the elements in one 
column made up of the elements in the previous column moved down one 
row. 


gust force. This equation represents a system of linear simultaneous 
equations where the order of the matrix is arbitrary; that is, the equa- 
tions may be written up to any desired value of s = m€. The solution 
for response can therefore be carried on as far as desired. Fortunately, 
the equations are of such a nature that simultaneous solution is not 

required. As mentioned, each of the matrices W, [E] , and [C] is 

triangular with all elements 0 above the main diagonal and with all 
elements on the main diagonal of each matrix equal; therefore, the main 

diagonal elements of [p] will also all have the same value and the 
elements above this diagonal will be 0. If each element on the main 


An equation in |p| alone is obtained by substituting |a | from 
this equation into equation (30) to yield 



(33) 


which is the basic response equation relating (3 (that is z-^") to the 


diagonal of [b] is denoted by d^_ and 

the main diagonal elements replaced by 0's, then 



D with 


[D] = d^l] + [ Dl ] 


With this equation, equation (33) may be written 



( 34 ) 


20 


NACA TN 2763 


Expanded, this equation has the form 



(35) 


It can be seen that a step-by-step solution for the successive values 

of 0 may now be made; that is, 0-l is solved for first then, with 0 L 

established, p 2 is solved for, and so on as far as is desired. With 

the value of |p| thus established, solution for |a| may now be made 
directly from equations (32). Values of the displacements z Q and z^ 
may be obtained directly from a and 0; z ± may be obtained from equa- 
tion (31); and z 0 may be obtained from this same equation with 0 
replaced by a. 

Some mention should be made with regard to the selection of the time 

interval 6 . A rough guide to use in selecting € can be obtained by 

considering X, which appears as the characteristic frequency in most 
response calculations. The period based on this frequency would be 

Pit 

T s = -j-. Experience has shown that a time interval in the neighborhood 

of 1/12 of this period yields very good results (in general less than 
1 percent error); accordingly, a reasonable guide in choosing e would 

be the equation e * — L ~. Some convenient value near that given by this 

C.K 

equation should be satisfactory; in general, it will be found that e 
may be 1 or greater. 

The procedure thus outlined provides a rather rapid evaluation of 
the response due to a prescribed gust. With the response thus evaluated 
the bending moment at any value of s or the complete time history of 
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bending moment may be found by application of equation (23). It should 
be evident that, if response values for either Zq" or ^l" sire known, 

the gust causing this response can be found by suitable manipulation of 
equations (30) and (32). Thus, if Zq" is known, (3 in equations (30b) 

and (32b) may be eliminated to give the equation 

{[A] * CotBlfc]' 1 ] l = | - |f| 


Direct substitution of z 0 " in this equation allows |f| to be deter- 
mined. In most practical cases the second term in equation (30b) con- 
tributes only a small amount and may be dropped with little resulting 
error in the gust force. The equation for |f| is then simply 

W 1=1 - |f| 


Determination of Response of One Airplane From 
Known Response of Another Airplane 

In general, a given gust condition produces different responses 
either for two different airplanes or for the same airplane with dif- 
ferent loading conditions or forward velocity. It would be expected, 
however, that the response of the two airplanes could be correlated 
through the common gust condition. This correlation may be demonstrated 
quite easily by means of the equations given in the preceding section. 
The case to be treated is as follows: The time history of bending 

moment due to a gust is assumed to have been measured in one airplane 
and it is desired to calculate directly from this time history what the 
bending moment due to the same gust would have been in another airplane. 
Although the derivation is presented in terms of bending moment, a 
similnr derivation could be made in terms of either accelerations or 
displacements . 

If use is made of equation (31) to write the successive values of 
the displacement z^ in terms of the accelerations, the bending-moment 

factor, equation (23), may be written in terms of the accelerations alone 
and the following matrix equation for K may be formulated: 
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Kl 


a l 


Pi 


k 2 


a 2 


p 2 


k 3 


a 3 


p 3 

0 

• 

= d 

• 

+ e 

• 

+ he 

• 


* 


• 




°m 


Pm 



1 

6 

1 

2 


1 

6 

1 i 


in - 1 


pi 

P 2 

p 3 




where 


4 . ,0 - ->0 

r l - r 2 


r l - r 3 

e = 2 - - m 

r l - r 2 


r-, - To o 

h = -±- 2- 

r l ~ r 2 


With the use of equations (32) this equation may be written 

Ik | = [e] I P I 


where 


[E] = 


— L c ] + e[l] + he 2 [g! 
^0 


(36a) 


(36b) 


(37) 
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in which Qf] is the identity matrix and 


[G] 


Substitution of the value of (3 as obtained from equations ( 36 ) into 
equation (33) gives the following relation between the gust forces and 
the bending-moment parameter: 



|f| = [D][E]“ 1 |K ffl | (38) 

The gust-force matrix |f| (see eq. (27b)) may be expressed in 
terms of the gust velocity by the following process: It is assumed that 

the initial vertical velocity of the gust is zero and that successive 
values of gust velocity of increment e are designated by u^, U 2 , 

u^ . • • • First-order difference equations are used to approximate the 
slope of the gust velocity, so that, in general. 


u , . w - V , i . 

m o a 


If this equation is used and the integral equation ( 27 b) is handled by 
the trapezoidal integration method similar to that used for equation ( 29 ), 
the gust force may be written in terms of the successive values of gust 
gradient so as to form the following matrix equation: 
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*1 

\[r 2 


*1 

^3 - 


^1 

t 2 





U 


u 2 


U 


u 3 


u 


u 4 


U 


(39) 


where i|^, \lr 2 , ^3 . . . are successive values of the \|f function. 

Substitution of this equation into equation ( 38 ) allows for the solution 
of — in terms of the parameter K as 


u 

U 


W' 1 WH' 1 KI 


(40) 


where [y] is the square matrix in equation (39)* Different airplanes 

flying through the same gust will experience the same vertical gust 
velocity for equal absolute distances of gust penetration; that is, 


(Vt) 


airplane 1 


= (Vt) 


airplane 2 


From equation (14), then, the following conditions must prevail: 


( SC o) 1 = ( sc o) 2 
( ec o) L = ( ec o) 2 


(41) 


where the subscripts 1 and 2 denote airplane 1 and 2, respectively. 
Satisfaction of the latter condition insures that the gust velocity as 
given by equation (40) would be the same for the two airplanes being 
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compared. This common gust condition may therefore he eliminated to 
yield the result 

(WTdIH^IkIX - (W- 1 [d][e]- 1 |k |) 2 


If it is assumed that K for airplane 1 is known, then K for air- 
plane 2 may he written 


K 


2 “ 


M [p] _1 m) 2 [[♦] _1 M M "Mki] 


(H£) 


where again the time interval chosen for the two airplanes must satisfy 
relation (4l). Thus, if the bending moment due to a given gust sequence 
is known for one airplane, the bending moment that would develop in 
another airplane encountering the same gust sequence can be determined 
from this known bending moment by the use of equation (42). If the mid- 
chords of the two airplanes are equal, the time interval may be taken 
equal and equation (42) reduces to 


|K| 2 =([E][D]- 1 ] 2 [MW- 1 |Ki] 1 (“3) 


SUMMARY OF CALCULATION PROCEDURE 


As a convenience, a summary of the basic steps necessary for calcu- 
lating the response of an airplane to a gust is given as follows: 

For accelerations and displacements: 


(1) With the use of the fundamental mode, wing plan form, and mass 
distribution, calculate the quantities Pq, p-^, X., r-^, and r ^ as 

given by equations (l 8 ). 


(2) Choose the time interval 



but for most cases 


e = 1 


6 . A convenient rule of thumb is 
should give satisfactory results. 


(3) Determine values of the unsteady-lift function 0=1-® at 
successive multiple intervals of e. (See fig. 1.) Also determine 
corresponding values of the gust-force integral f(s), equation ( 27 b). 

As an aid, curves for f(s) are presented in figure 1 for the sharp-edge 
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p 

gust and in figure 2 for various-length sine gusts, sine gusts, and 
triangular gusts. (The curves in fig. 1 have been obtained from 
eqs. (9) and (10). These approximations, although rather accurate for 
the lower values of s, are noted to cross; actually, they should not 

cross and are known to have the same asymptotic approach to unity. ) 

(4) From the following definitions: 


Af - d 0 + e0 O 


Am = 2e0 m-l 


(m > 1) 


Bn = 


r i e0 o 


Sm = 


2r 60 
1 m-1 


: i-^r( 1 + 4 £ ) + (^- r i) ee o 


(m > 1) 




- (« - 1)^ A 2 - 2 - rd c Vl 


(m > 1) 


set up the following matrices: 


[Aj = 


A-2 A 1 


A 3 A 2 A 1 


A^ A^ Ag A^ 
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[B] = 


Bn 


B 2 B 1 


B3 Bg B^ 


Bi, Bo Bo B , 


[C] - 


C 2 C 1 


C^_ 

c 4 c 3 c 2 C 1 


Then calculate the matrix 


[D] - 


^[A][C] ♦ [B] 
^0 


(5) Solve for the values of P ( which equals z ]_") from equa- 
tion (33), by the method outlined following equation (33)- (See 

eq. (3*0* ) The values of Zj_ and a ^ which equals Zq" ) can then be 

calculated from equations (31) and (3^). 

For bending moment: 

(6) In order to compute bending moment, determine r^, t}q, and t^ 

as given by equations (24), where M^, M c ^, anc ^ ^ C 1 in ^ ese 

equations depend on the particular wing station being considered and are 
given by equations (21). 
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(7) Determine bending moment by use of equation ( 23 ) with the 
values of response already established. This equation may be applied 
directly to any desired time value. Maximum bending moment usually 
occurs very close to the time when is a maximum. 

EXAMPLES 


Example A .- In order to provide an illustration and give an idea 
of the accuracy of the present analysis, the response to a sharp-edge 
gust of the two-engine-airplane example considered in reference 1 was 
determined. The weight distribution over the semispan, the wing-chord 
distribution, and the fundamental bending mode are shown in figures 3 , 

4, and 5 . The frequency and deflection of the fundamental mode were 
calculated by the method given in reference 10. The solution is made 
for a forward velocity of 210 miles per hour and a gust velocity of 
10 feet per second. 

The lift-curve slope used in reference 1 was 5.4l; to be consistent, 
the same value was used here. Furthermore, the unsteady-lift function 
used for a change in angle of attack in the example presented in refer- 
ence 1 was given by the equation 


(1 - ®) a=6 = 1 - 0.36le-°*38ls 

rather than by equation ( 9 ). Thus, this equation was also used here. 

The gust unsteady-lift function used was that given by equation (10). 

The various physical constants and the basic response and bending- 
moment parameters are given in table 1 ; the values of the unsteady-lift 
function and the values of the gust force are listed in table 2. The 

matrices [A], M , and [c] used in the solution are given in table 3 • 

The solution for response is shown in figure 6 (a) where the deflec- 
tion coefficients aQ and a-^ in inches are plotted against distance 

traveled in half-chords. The corresponding deflection quantities for 
the example given in reference 1 were determined and, for comparison, 
are also shown in the figure. A similar comparison is made in fig- 
ure 6 (b) for bending stresses at the fuselage and engine stations, 
stations 0 and 1 from reference 1. The agreement is seen to be good. 

Example B .- A second example is included in order to illustrate one 
means by which the method may be used to evaluate the influence of bending 
flexibility upon the response to a gust. The physical characteristics for 
the airplane considered in this example are listed also in table 1 , and 
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equation (9) is used for the function 1 - $ instead of the values 
given in table 2. Maximum values of the bending moment that develops 
at the fuselage station during flights through sine gusts of various 
lengths have been determined, both for the airplane considered flexible 
and for the airplane considered rigid. The results are shown in fig- 
ure 7(9*) where maximum values of* "the bending-momeirt parameter K are 
plotted against gust-gradient distance H. The difference between the 
two curves represents the increase in bending moment due to effects of 
wing bending flexibility. By taking the ratio of K for the flexible 
case to K for the rigid case, a type of dynamic response factor is 
formulated which gives a direct measure of the influence of wing flexi- 
bility. This ratio is designated 7 M and is shown in figure 7(b). As 

an example of the significance of this plot, the value of / M = 1.16 

® 7 ^ means that flexibility results in a 16— percent dynamic over- 
shoot in the stress from the value that would be obtained at H = 5 on 
the basis of a rigid-body analysis. It may be seen also that the value 
of r M is approximately unity for values of H = 10 and greater; there- 
fore, in this range of gust-gradient distances a rigid-body treatment 
would be sufficient for this airplane. 


DISCUSSION 


The derivation presented herein is intended to provide a convenient 
engineering method for calculating the response of an airplane to a gust 
where wing bending flexibility is included. The method is believed to be 
well-suited for making trend studies which evaluate, for example, the 
effect on response of such factors as mass distribution, speed, and 
altitude. Although the unsteady-lift functions for two-dimensional 
unsteady flow are presented, the method is general enough so that the 
unsteady-lift functions for finite aspect ratio, for subsonic compress- 
ible flow, and for supersonic flow may be used as well. (See refs 7 
and 11 to 15. ) v . 1 

Since the numerical method is based on an integration procedure, it 
possesses the desirable feature that a fairly large time interval may be 
used and good accuracy still obtained. As an accuracy test, solutions of 
equations (16) and (17) were made for several cases by the exact Laplace 
transform method as well as by the numerical process, in which process 
the time interval was selected according to the rule of thumb suggested. 
When the results were plotted to three figures, the difference between 
the two solutions was barely discernible. 

Additional bending modes could be included in the analysis but this 
refinement is really not warranted. Some calculations made with addi- 
tional modes gave results which differed only slightly from the results 
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obtained when only the fundamental mode was used. The good agreement 
of results obtained in example A with the results obtained by the more 
precise method given in reference 1 also illustrates this point. 
Furthermore, if additional degrees of freedom are to be used, it would 
appear more important to include wing torsion and airplane pitch. The 
extent to which torsion influences the results is probably governed most 
by the nearness to the flutter speed. The importance of airplane pitch 
is probably governed most by the gust length} some investigations dealing 
with pitch have indicated that except for very light wing loadings the 
pitch of the airplane does not influence the results appreciably until 
gust lengths of from 20 to 30 chords or larger are involved. Thus, the 
present analysis, although limited to the degrees of freedom of vertical 
motion and wing bending, should probably be sufficiently satisfactory 
for speeds near the cruising speed and for gust-gradient distances up to 
approximately 10 chords. 

The analysis may be useful in assessing the significance of wing 
flexibility in the technique of measuring gust intensity by means of an 
airplane. In this technique gust severity is usually measured by means 
of an accelerometer placed at the center line of the airplane. In order 
to obtain a rough idea of whether flexibility may have some effect on 
this measurement, calculations for the maximum accelerations at the 
center line and for the maximum acceleration at the nodal points (the 
true center- of- gravity acceleration) may be made for various assumed 
gust lengths. A comparison of these computed maximum acceleration values 
should give some idea as to the extent to which wing flexibility may 
alter the measurements in actual flight. 


CONCLUDING REMARKS 


The analysis presented herein for the response of an airplane to a 
gust should provide a useful means for evaluating the effects of wing 
flexibility. A convenient, but accurate, numerical solution of the 
response equations is developed which is well-suited for trend studies 
such as the evaluation of the effects of mass distribution, speed, 
altitude, and similar factors. 

As indicated by an example, the method gives good agreement with 
the results of the more precise but more lengthy recurrence matrix 
method of NACA Rep. 1010. 

The method permits the evaluation of a gust causing a known response. 
A procedure is given wherein the known gust response of one airplane may 
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be used directly to determine what the response would be for another 
airplane flying through the same gust condition. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 21, 1952. 
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APPENDIX 


DERIVATION OF EQUATION RELATING DISPLACEMENT TO 
PREVIOUS SUCCESSIVE VALUES OF ACCELERATION 


In this appendix, a derivation is given of equation (31) which gives 
the value of displacement in terms of successive past-history values of 
acceleration. Suppose that the second derivative (acceleration) of a 
function is approximated by a succession of straight-line segments as 
shown in the following sketch: 



where the segments cover equal 
initial condition that Zq = 0 


intervals € of the abscissa s and the 
is assumed to apply. If a dummy origin 


is now considered at the station m - 1, the segment between sta- 
tions m - 1 and m may be represented by the equation 


z 


11 


11 


z 


m-1 


+ 




1 


€ 


S 


Two successive integrations give the relations for z* m and z m as 
follows: 


z 1 


in 


- z 


m-1 ^2 


m-r 


2 € 


z 


t 

m-1 


q 2 

_ 1 » o , 

z ~ z m-1 p 


m 


- z 


m-1 


S 3 + 


m-1 


s + z 


m-1 


5E 
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where the constants of integration z, m_i an d (initial conditions 

for the interval) have been introduced. If s is set equal to e 
these two equations, the following equations result: 


in 


m 


> !(*"„ ♦ z "m-l) 


+ z 


m-1 


(Al) 


Z™ = ~ 

m £ 


m 


m-1 


+ z' 


m-l e 


+ z. 


’m-1 


(A2) 


From these two equations the values of z’ m and z m at any time 

interval may be given in terms of the second derivative at all previous 
time intervals. For example, with initial conditions of z q = z’q = 0, 

equation (Al) becomes for m = 1 


z 


1 



(A3) 


and. for m = 2 


z 


2 




z 


1 


Combining this equation and. equation (A3) results in the relation 


z ’ 2 = e ( z "l + \ z " 2 ) 

This process may be carried, through for each of the time stations to 
yield, the following general equation for z* m : 


z 


1 

m 


efz\ + z" 2 + z" 3 + 


+ Z 


m-1 


1 z ” \ 

2 Z my 


(A4) 


which, of course, is the trapezoidal approximation of the area under the 
z"-curve. Equation (A2) for z m may be treated similarly and it is 
found that the general equation for Zju may be written 


„2 


J m 


(m - l)z"^ + (m - 2 )z"2 + . 


2z m-2 + 2 


> 1 " 
m-1 + ^ 2 m 


(A5) 
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This equation thus gives the displacement at any time station in terms 
of the accelerations at all previous time stations. 

It may he noted that, if higher-order segments (parabolic or cubic) 
had been used instead of straight-line segments to approximate the 
second derivative, equations similar in form to equations (A4) and (A5) 
would also result. For most practical purposes, however, the accuracy 
of equation (A5) is sufficiently good as long as the interval € is 
chosen so that the straight-line segments roughly approximate the second 
derivative. 
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TABLE 1.- PHYSICAL CHARACTERISTICS FOR EXAMPLE AIRPLANES 


Example A 


Example B 


37,450 

870 

1120 

154 

0.0765 

308 

10 

0.0208 

1.0 

5.41 

64.16 

0.9045 

0.4353 
0. 2181 

0.1358 

0.452 

0.547 

23.49 

10.19 

3.665 

3.391 

0 . 00537 
0.00669 


33,450 

870 

1120 

164 

0.0765 

374 

0.0183 

1.0 

6.28 

46.8 

0.748 

0.392 

0.225 

0.143 

0.457 


15.94 


2.555 


0.00537 


w, lb 
s, ft 2 

b, in. 


'O’ 


in. 


, lb/ ft 3 . . . 
V, ft/sec . . . 
U, ft/sec . . . 

Jsec . . . . 
[half- chords 


K 0 

^1 
\ . 

r l 


r 3 


Jfuselage station 
[engine station 

« Jfuselage station 
(engine station 

Jfuselage station 
^engine station 


= , in 


-3 jfuselage station 
[engine station 


here denotes distance from neutral axis to extreme fiber. 
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TABLE 2.- 1 - <t> ORDINATES AND GUST-FORCE ORDINATES 

FOR SHARP-EDGE GUST, 6 = 1.0 


m 

9 m or (1 - <k ) a=6 

f or ^ 

0 

0.6390 

0 

1 

.7534 

.377 

2 

.8315 

.547 

3 

.8849 

.635 

4 

.9214 

.692 

5 

.9463 

.736 

6 

.9633 

.771 

7 

.9749 

• 798 

8 

.9829 

.821 

9 

.9883 

.845 
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TABLE 3.- MATRICES USED IK EXAMPLE A 


A Matrix 


64. 799 

1.5068 61*. 799 

1.6630 1.5068 

1.7698 1.6630 
1.81*28 1.7698 
1.8926 1.81*28 
1.9266 1.8926 
1.9498 1.9266 
1.9658 1.91*98 
1.9766 1.9658 


64.799 

1.5068 64. 799 

1.6630 1.5068 

1.7698 1.6630 
1.81*28 1.7698 
1.8926 1.8428 
1.9266 1.8926 
1.9498 1.9266 


64.799 

1.5068 64.799 

1.6630 1.5068 

1.7698 1.6630 
1.8428 1.7698 
1.8926 1.8428 


64. 799 

1.5068 64.799 
1.6630 1.5068 
1.7698 1.6630 


64 . 799 

1.5068 64.799 


B Matrix 


0.1394 

.3286 O.I394 

.3627 .3286 

. 3860 . 3627 

. 4019 . 3860 

. 4128 . 4019 

. 4202 . 4128 

. 4252 . 4202 

. 4287 . 4252 

.4311 .4287 


0. 1394 

.3286 0.1394 
.3627 .3286 
. 3860 . 3627 
.4019 .3860 
. 4128 . 4019 
. 4202 . 4128 
. 4252 . 4202 


0.1394 


.3286 

0.1394 

.3627 

.3286 

.3860 

.3627 

.4019 

.3860 

.4128 

.4019 


0.1394 

.3286 0.1394 
.3627 .3286 
. 3860 . 3627 


0.1394 
• 3286 


0.1394 


C Matrix 


4.5367 

1.3954 4.5367 

2.2445 1-3954 
3.0735 2.2445 

3.8889 3.0735 

4.6949 3.8889 
5. 4947 4. 6949 
6.2900 5.4947 
7.0824 6.2900 
7.8726 7.0824 


4.5367 

1-3954 4.5367 

2.2445 1.3954 
3.0735 2.2445 

3.8889 3.0735 

4.6949 3.8889 
5.4947 4.6949 
6. 2900 5. 4947 


4.5367 

1-3954 4.5367 

2.2445 1.3954 
3-0735 2.2445 

3.8889 3.0735 

4.6949 3.8889 


4.5367 

1.3954 4.5367 
2.2445 1.3954 
3.0735 2.2445 


4.5367 

1.3954 4.5367 



Figure I r Unsteady- lift functions (see equations (9) and (10) ) where, for a sharp-edge gust, 

the gust force f(s)=\|/(s). 
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Sin gust 
Sin 2 gust 
Triangle gust 



(a) H=2.5 chords. 

Figure 2.- Value of the gust-force integral f (s ) = J -g- V(s-cr)dcr 


for three gust shapes. 
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Figure 2 Continued. 
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(c)H= 7.5 chords. 
Figure 2.- Continued. 
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s, Half-chords 
(d)H=IO chords. 
Figure 2 Concluded. 
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Figure 3r Semispan weight distribution for the two-engine airplane of 

example A. 



Figure 4.- Wing chord distribution for airplane of example A. 
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Present analysis 

Reference I 


a 0 and a,, 
in. 



Stress, 



(b) Stresses. 


Figure 6 r Response of example airplane A to a 10-foot-per-second sharp-edge gust. 

V=2IO miles per hour. 
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Figure 7 r Bending moment and dynamic response factor for airplane 
of example B due to flight through sine gusts. 
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